Manifestation of Topological Protection in Transport Properties of Epitaxial Bi2Se3 

Thin Films 
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The massless Dirac fermions residing on the surface of three-dimensional topological insulators 
are protected from backscattering and cannot be localized by disorder, but such protection can be 
lifted in ultrathin films when the three-dimensionality is lost. By measuring the Shubnikov-de Haas 
oscillations in a series of high-quality Bi2Se3 thin films, we revealed a systematic evolution of the 
surface conductance as a function of thickness and found a striking manifestation of the topological 
protection: The metallic surface transport abruptly diminishes below the critical thickness of ~6 
nm, at which an energy gap opens in the surface state and the Dirac fermions become massive. At 
the same time, the weak antilocalization behavior is found to weaken in the gapped phase due to 
the loss of TT Berry phase. 

PACS numbers: 73.25.-|-i, 71.18.+y, 73.20.At, 72.20.My 
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In topological insulators (TIs) the energy states are 
fundamentally modified from ordinary insulators by 
strong spin-orbit interactions, giving rise to a topolog- 
ically distinct state of matter with a gapped insulating 
bulk and a gapless metallic surface pLj]. Various inter- 
esting phenomena, including surface transport of spin- 
filtered Dirac fermions that are immune to localization, 
have been predicted and raised expectations for novel 
applications 0-01 ■ However, the progress in real applica- 
tions of TIs crucially relies on the ability to manipulate 
the surface current in transport experiments. At present, 
such basic characterization as the surface conductance 
measurement has been possible only in a few cases in 
single crystals [sl-fl^ because of the dominance of bulk 
transport caused by unintentional doping due to defects. 
Molecular beam epitaxy (MBE) is a promising technique 
for the synthesis of TIs |13l4lg | owing, in part, to the rela- 
tively low deposition temperature at which defect concen- 
trations can be reduced from those in bulk crystals grown 
in thermal equilibrium. So far, using MBE-grown films, 
angle-resolved photoemission spectroscopy (ARPES) (loj 
and scanning tunneling spectroscopy (STS) have pro- 
vided useful information about the topological surface 
state (SS), and in transport experiments, such phenom- 
ena as weak antilocalization (WAL) and gate -controlled 



ambipolar transport have been reported [21|, |24|. Also, 



since the MBE technique gives a precise control over film 
thickness, transport measurements for widely varying 
surface-to-bulk conductivity ratio have been performed 
13- 25 1, although a reliable separation of surface Dirac 
electrons from bulk carriers has been hindered by a rel- 
atively low mobility of carriers in available thin films. 
Recently, we have succeeded in growing high-quality epi- 
taxial films of Bi2Se3 that have a sufficiently high sur- 
face electron mobility to present pronounced Shubnikov- 
de Haas (SdH) oscillations. This made it possible to di- 
rectly probe the surface conductance and the topological 
protection of the SS. 



The immunity of the surface Dirac fermions to localiza- 
tion has a twofold origin [IHl] ■ One is the it Berry phase 
associated with massless Dirac fermions [2^, which pro- 
tects them from weak localization through destructive in- 
terference of time-reversed paths. The other is the pecu- 
liar spin-momentum locking which nulls the backscatter- 
ing probability [H-Q- Those mechanisms are collectively 
called topological protection. Recently, it was found 
(lol I27l - l3l| that when TIs are thinned to the extent that 
the top and bottom surface states feel each other, their 
hybridization leads to an opening of the gap at the Dirac 
point and results in a degenerate, massive Dirac disper- 
sion. This gapped state obviously violates the topological 
protection, but its consequence in the surface transport 
properties has not been duly addressed. In fact, this 
question is important because recently a lot of attention 
has been paid to the way to open a gap at the Dirac 
point 32| to realize topological magnetoelectric effects 



[33|. One may expect that unless the chemical potential 
is located exactly within the gap, the metallic surface 
transport is largely unaffected by the gap opening be- 
cause states are kept being available for transport at the 
Fermi level. However, in the present work it turned out 
that the change in the Dirac spectrum deep in the oc- 
cupied state has a profound effect on the physics at the 
Fermi level. 

The growth of Bi2Se3 films occurs in a layer- by- layer 
manner, in which the 0.95-nm thick Se-Bi-Se-Bi-Se quin- 
tuple layer (QL) constitutes the basic unit jl^. Our 
systematic magnetotransport measurements for varying 
thickness reveal a sudden diminishment of the surface 
transport below the critical thickness of ~6 QL, below 
which the energy gap opens in the Dirac spectrum jl9| . 
We also observed that the weak antilocalization behav- 
quickly weakens below the critical thickness. 
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We discuss that those striking effects are due to acquired 
degeneracy of the surface states [2^ and loss of their tt 
Berry phase [34|, |35| in the gapped phase. 
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FIG. 1: (Color online) (a) AFM image of a 50-nm thick Bi2Se3 
film showing atomically flat terraces with 1-QL steps. Inset: 
typical RHEED pattern, (b) X-ray diffraction pattern of a 
200-QL film, (c) Temperature dependences of Rs for different 
thickness. 



Our MBE films were grown under Se-rich conditions on 
insulating sapphire (0001) substrates whose size was ap- 
proximately 15 X 4 mm^. To obtain films of high enough 
quality to present SdH oscillations, we employed a two- 
step deposition procedure [13, Ell • Both Bi (99.9999%) 
and Se (99.999%) were evaporated from standard Knud- 
sen cells. The Se2(Se4)/Bi flux ratio was kept between 
15 - 20. The growth rate, which is determined by the Bi 
flux, was kept at 0.2 - 0.3 QL/min. The resistance R^x 
and the Hall resistance Ryx were measured in the Hall 



bar geometry by a standard six-probe method on rect- 
angular samples on which the contacts were made with 
silver paste at the perimeter and cured at room tem- 
perature under pressure of '^1 Pa. The magnetic field 
was swept between ± 14 T at fixed temperatures and 
was always applied perpendicular to the films, except for 
the angular-dependence measurements of the SdH oscil- 
lations. 

An atomic force microscopy (AFM) image of our rela- 
tively thick (50 QL) film is shown in Fig. 1(a), where 
a large, atomically fiat area of ^1 /im^ and several 
sharp terraces with the height of exactly 1 QL can be 
clearly seen. The reflection high-energy electron diffrac- 
tion (RHEED) pattern with sharp 1x1 streaks [Fig. 1(a) 
inset] and sharp x-ray diffraction peaks [Fig. 1(b)] are 
also indications of the high crystal quality of our films. 

Temperature dependencies of the sheet resistance, 
Rs{T), measured in films with systematically changed 
thickness t down to 2 QL are shown in Fig. 1(c). The 
Rs{T) behavior is metallic in thick films, but below t = 5 
QL, it starts to show an upturn at low temperatures. In 
particular, the sharp divergence in Rs for T ^ in the 
2-QL film is indicative of strong Anderson localization 
and an insulating ground state [37 1. 

The breakthrough in the present work is that our films 
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FIG. 2: (Color online) Surface SdH oscillations in the 10-QL 
film, (a) dRyx/dB in tilted magnetic fields, plotted as a func- 
tion of 1/B± {— 1/B cos6); curves are shifted vertically for 
clarity. Dashed lines mark the positions of maxima. Inset 
shows the geometry of the experiment, (b) Dingle plot of 
the oscillations in ARxx at 1.6 K, obtained after subtracting 
a smooth background from Rxx{B), giving the Dingle tem- 
perature of 16 K. Inset: T-dependence of the SdH amplitude 
for 6 = 0°. (c) Landau- level fan diagram for oscillations in 
Gxx measured at T = 1.6 K and e = 0°; following Ref. 
integers n (half- integers n + i) are assigned to the minima 
(maxima) in AGxx- The solid line is a linear fitting to the 
data with the slope fixed at F = 106.8 T; the dashed line has 
the same slope and extrapolates to zero. Upper inset shows 
AGxx vs. 1/B after subtracting a smooth background; lower 
inset shows its Fourier transform giving F — 106.8 T. (d) 
Fitting of the two-band model to the Ryx{B) data at 1.6 K. 



exhibit pronounced two-dimensional (2D) SdH oscilla- 
tions to provide a direct way to probe the surface charge 
transport. As an example, the analysis of the SdH os- 
cillations in the 10-QL film is shown in Fig. 2. The 2D 
character of the oscillations is evident in Fig. 2(a), where 
the positions of the maxima and minima depend only on 
the perpendicular component of the magnetic field, B±. 
The oscillation frequency F = 106.8 T is obtained from 
the Fourier transform [lower inset of Fig. 2(c)], and this 
is a direct measure of the Fermi wave number kp = 5.7 
X 10*^ cm" 
tal Material 



^ As we discuss in detail in the Supplemen- 
37|, if the SdH oscillations are due to the 



trivial 2D electron gas (2DEG) which may form due to 
a band bending near the surface '38^, this kp is so small 
that it imposes too strong a constraint on the possible 
bulk Fermi level, which makes it impossible to consis- 
tently explain the transport data. Hence, we identify the 
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oscillations to be due to surface Dirac fermions, and the 
obtained kp gives their density Us = 2.6 x 10^^ cm~^. 

In this 10-QL film, we observed only a single frequency, 
but we usually see two frequencies in other films (see 
Fig. S3 of the Supplemental Material [ll]), suggesting 
that the top and bottom surfaces have somewhat differ- 
ent Us- The temperature dependence of the oscillation 
amplitude [Fig. 2(b) inset] gives the cyclotron mass rric 
= 0.2mo (mo is the free electron mass) [39| which in turn 
gives the Fermi velocity vp = 3.3 x lO'' cms~^. This 
vp is consistent with the ARPES data [i^ as well as the 
STS data ^ Slj for the Dirac cone. The obtained kp 
value corresponds to the Fermi level of ~0.16 eV above 
the Dirac point for the topological surface state, which 
points to a slight upward band bending 42| . The Dingle 
analysis [Fig. 2(b), see the Supplemental Material [37| 
for details] yields the mobility = 1330 cm^V~^s~-^. 
Finally, Fig. 2(c) shows the Landau-level fan diagram 
for the oscillations in conductance Gxx, where the posi- 
tions of the minima in /S.Gxx (shown in the upper inset) 
are plotted as a function of n [431 ]. Here, to minimize the 
error occurring from extrapolation [l^ . we fix the slope 
of the linear fitting by using F = 106.8 T obtained from 
the Fourier analysis and determine the intercept n = 0.40 
±0.04 (solid line); this is very close to the ideal value of 
0.5 for Dirac electrons bearing the tt Berry phase [3], 
giving further confidence in the origin of the SdH oscil- 
lations. For comparison, a straight line with the same 
slope to give zero Berry phase is shown in Fig. 2(c) with 
a dashed line, which is obviously inconsistent with the 
experimental data. 

To estimate the contribution of the SS in the overall 
transport in this 10-QL film, we use the magnetic-field 
dependence of the Hall resistivity, Ryx{B) [Fig. 2(d)], 
which is not linear in B and thus signifies the presence 
of at least two types of carriers. The fitting of a stan- 
dard two-band model [^, 0, 11^ to the data, in which 
is fixed by the SdH frequency, gives the surface contribu- 
tion to the total conductance, Gs/Gtot^ of 28%. The /i^ 
value obtained from this fitting is close to the SdH result, 
assuring the consistency of our analysis. 

The same analysis can be applied to all measured films 
with t > % QL, in which we consistently observed SdH 
oscillations [37|. Evolutions of the transport parameters 
with changing t are summarized in Figs. 3(a)-3(c). We 
note that by tracing the evolution of the SdH oscillations 
starting from thick films, we can distinguish the topo- 
logical SS from the 2D quantum-well (QW) state of the 
bulk origin [37| as the source of the SdH oscillations. 

Our main finding is that the surface transport abruptly 
diminishes below a critical thickness tc which is located 
between 5 and 8 QL. This change is most convincingly 
manifested in the behavior of Ryx{B), which suddenly 
becomes i3-linear in films with t < 5 QL [Fig. 3(d)]; 
this indicates that the transport becomes suddenly dom- 
inated by only one type of carriers. Correspondingly, 
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FIG. 3: (Color online) (a) Squares are the Us obtained from 
SdH oscillations (open and filled squares represent two dif- 
ferent surfaces of the same film) and the dashed squares are 
an extrapolation of the trend above ta\ circles are the sheet 
density of bulk carriers obtained from two-band analyses, (b) 
Mobilities of surface (squares) and bulk (circles) carriers ob- 
tained from two-band analyses, (c) Gs/Gtot obtained from 
the data in (a) and (b). Insets are schematic pictures of the 
energy bands for the two regimes, (d) Ryx{B) at 1.6 K for 
various thickness shown in QL unit; inset shows the derivative 
dRyx/dB for the 5 and 8 QL films (5-QL data are shifted for 
clarity), (e) Rxx[B) / Rxx[Q) of the same films. 



the SdH oscillations disappear for t < tc- More quantita- 
tively, assuming that is essentially unchanged through 
tc [dashed squares in Fig. 3(a)], one can estimate that 
fis must be suddenly degraded by more than an order 
of magnitude in samples with t < tc [dashed squares in 
Fig. 3(b)] for R„x(B) to become linear and be governed 
by bulk carriers |37| . The magnetoresistance behavior 
shown in Fig. 3(e) also presents a qualitative change be- 
low tc, showing a negative slope at high fields. This evo- 
lution is best represented in the t dependence of Gs/Gtot 
[Fig. 3(c)], which shows a steady increase with decreas- 
ing t to reflect the change in the surface-to-bulk ratio, 
but it drops sharply below tc to signify that the surface 
transport is abruptly diminished. 

This observation naturally calls for the question 
whether the observed diminishment of the surface trans- 
port in ultrathin films might be related to a lowering of 
the quality in thinner films. In this respect, our ultrathin 
films remain essentially flat and smooth across tc, with 
the surface bumpiness of only '^1 QL (see Fig. SIO of 
the Supplemental Material [37|). This observation, com- 
bined with the fact that Gs/Gtot increases steadily with 
decreasing t until it reaches tc, testifies against the above 
concern. 

In addition to the above results, we found a striking 
change in the WAL behavior [21, 22| below t = 5 QL. Fig- 
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FIG. 4: (Color online) (a) WAL behavior in sheet conductance 
at 1.6 K for various thicknesses shown in QL units; dashed 
lines are the fittings using Eq. (1). (b) Thickness dependence 
of a. Inset shows schematic energy bands above and below 
the critical thickness. 



ure 4(a) shows the magnetoconductance of our films mea- 
sured in perpendicular magnetic fields at 1.6 K. Dashed 
hues are the fitting with the Hikami-Larkin-Nagaoka for- 
mula 
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AC(B) =a 



In 



4eL^B 2 / y^eL^B I 

(1) 

where ^' is the digamma function and is the phase 
coherence length. The prefactor a should be — ^ for each 
transport channel that either carries a tt Berry phase 
p6| or bears a strong spin-orbit interaction [4^. In our 
analysis, a and are the only fitting parameters and 
Fig. 4(b) shows the t dependence of obtained a (see the 
Supplemental Material [37| for details). For t > 5 QL, 
we observed a « — ^ similar to that reported for metal- 
lic Bi2Se3 thin films where top and bottom surfaces are 
connected though bulk electrons 2l|, |22[ . The change to 
a w observed for t < 3 QL is in accord with the dimin- 
ishment of the surface transport channel and the even- 
tual localization of the bulk state. A similar tendency 
was also observed in previous studies 23-25|, although 
the decrease in a was less pronounced probably due to a 
larger metallicity of the measured samples. 

Now we discuss the origin of our observations. As was 
already found in Bi2Se3 ultrathin films by ARPES To*], 
an energy gap in the SS opens at the Dirac point below 
tc ~ 6 QL. Obviously, our transport measurements reflect 
this change in the Dirac dispersion. The gap is due to hy- 
bridization between top and bottom surfaces [schemati- 
cally shown in the Fig. 4(b) inset], and such a hybridiza- 
tion gap A changes the massless Dirac dispersion E — 
ihvpk into a massive one E = ±\/ {hvpkY + (A/2)-^ 
Ulail. For BizSea, A is about 0.25 eV at t = 2 QL 
In our films, the Fermi level is estimated to be ~0.16 
eV above the Dirac point, so even in our thinnest film 



the Fermi level is not in the gap but crosses the SS as 
schematically shown in the inset of Fig. 4(b). Hence, 
the observed drastic suppression of the surface transport 
is not due to the disappearance of surface carriers but is 
likely due to an enhanced scattering of the carriers [i^ . 

The gap opening also has a profound effect on the 
Berry phase 0s of the surface band ipkir)- In the sim- 
plest case, (/)_B is given by [34 1 



dip ( ipkir) 



TT 1 



and hence the Berry phase is reduced from tt when a 
gap opens. This lifts the immunity of the SS from weak 
localization. Also, this change in 0b weakens the WAL in 
the SS [35]. It is useful to note that in Fig. 4(a) the WAL 
was still observed for i = 5 QL which is below ic, but this 
is natural because the small A at this t [l^ makes the 
deviation of 0b from tt to be small. 

More importantly, in the gapped phase, the SS be- 
comes degenerate |29|], which means that now for each 
momentum both up and down spin states are available. 
This opens the backscattering channel and significantly 
reduces the surface mobility. Actually, the hybridiza- 
tion of top and bottom surfaces means that the system 
is no longer truly three-dimensional (3D), so it is natu- 
ral that the topological properties of 3D TIs are lost in 
the gapped phase, in full agreement with recent ah ini- 
tio density functional studies of TI thin films [s^, Isij . 
In this respect, the present observation is a spectacular 
manifestation of the topological protection of the SS in 
3D TIs. 
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Supplemental Material 



SI. Van der Waals epitaxy and two-step deposition 

The building block of the layered 612863 is a quintuple layer (QL) as depicted in Fig. Sl(a). The surface is 
terminated by Se atoms without dangling bonds, interacting only by weak van der Waals forces. Due to the weak 
interaction with a substrate, an epitaxial growth with its own lattice constant is possible at the very beginning (with 
a drastically relaxed lattice matching condition) as schematically shown in Fig. Sl(b). Thus, the growth of ultrathin 
films starting from one QL is possible [48]. This is called van der Waals epitaxy. 

One of the key parameters of the MBE growth is the temperature of deposition, which determines the rate of ad- 
sorption/desorption of Bi and Se as well as their surface diffusion. We found that increasing the substrate temperature 
greatly improves the size of grown terraces. Unfortunately, at high temperatures, the growth of thin films is virtually 
impossible due to a weak interaction of adatoms with the substrate. The crucial improvement was achieved by using 
a two-step deposition method [i^: In the first step, the initial 1 QL of Bi2Se3 is deposited at a low temperature of 
110-120°C; then, the temperature of the substrate is slowly raised up to 300-320°C, and, in the second step, the 
desired number of QLs is deposited. We obtained the best quality thin films at the temperature of the second step of 
around 300°C. 



S2. Anderson localization 

An upturn in the sheet resistance in ultrathin Bi2Se3 films has always been seen at low temperatures [ii-il, m 
spite of the difference in concentration of charge carriers and their mobilities. Figure S2 shows the diverging trend of 
Rs at low temperatures we observed in the 2-QL film. A rapid decrease in the mobility with decreasing temperature 
is also shown in the inset of Fig. S2, which points to a localization of charge carriers in the zero-temperature limit. 



Similar behavior of RsiT) was studied in Ref. 2j| on films with a higher conductance, where the insulating ground 
state was proposed to be due to strong interactions of surface electrons, which were assumed to dominate the transport 
in ultrathin films. However, for ultrathin films with lower conductance, which is determined by the concentration 
of charge carriers and the strength of disorder, the Anderson localization of a disordered 2D metallic conductor 
seems to be a better (more natural) explanation for the observed insulating state. Such localization behavior in low- 
conductance ultrathin films highlights the disappearance of topological protection of the surface state. This conclusion 
is also supported by a recent study of field-effect devices fabricated from ultrathin (^3.5 nm) Bi2Se3 crystals obtained 
by mechanical exfoliation on Si02/Si substrates [s^l- By measuring the temperature and gate- voltage dependences 
of the conductance, the authors of Ref. [H^l observed clear insulating behavior with an activated energy gap when 
the chemical potential is tuned into the gap, which made them conclude that their 3.5-nm thick Bi2Se3 crystals are 
conventional insulators. 



S3. SdH oscillations 



Transport properties of surface carriers are obtained mostly form the SdH oscillations observed in films with the 
thickness t above 6 QL (Fig. S3). The frequency F identified in SdH oscillations in i? vs. 1/B gives the size of a 2D 
Fermi surface (FS) via the Onsager relation F — {h/2Tre)A, where A — irkp is the FS cross section; h and e are the 
reduced Plank's constant and the electron charge, respectively; kp is the Fermi wave number. The concentration of 
spin-filtered surface charge carriers is — kp/Air. 

The cyclotron mass rUc is obtained from the temperature damping factor of the first harmonic, Rt — ttA/ sinh(7rA), 
where A — 2i:kBT /hujc and Wc = eB/rric is the cyclotron frequency From kp and mc, we obtain the Fermi 

velocity vp = hkp/nic. 

The quantum scattering time r = h/{2TTkBTp)) and the mean free path is = vpT are obtained from the Dingle 
damping factor Ro = exp(— 27r^A:B?D/?ia;c), which accounts for scattering effects [39| : the surface mobility /is is 
estimated from the relation for the surface conductance Gg — efigris — {e'^ /h){kp£s). 

S4. Sudden drop in the surface mobility below tc 

As has been discussed in the main text, the SdH oscillations disappear below which coincides with the disap- 
pearance of topological protection of the surface state (SS) when the energy gap opens and the Dirac fermions become 
massive. Here we present additional discussions that the disappearance of the SdH oscillations with decreasing film 
thickness cannot be simply due to a smooth decrease of the surface mobility below the observational limit. 
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Figures S4 and S5 show Rs{B) and Ryx{B) data measured in films with t = 9> and 5 QL, respectively, together with 
their fittings considering the surface and bulk channels, similar to what has been described in the main text for t = 
10 QL [Fig. 2(d)]. For the 8-QL film (Fig. S4), where we observed two SdH frequencies (Fig. S3), the concentrations 
of surface carriers ni = 1.7 x 10^^ cm~^ and n2 = 5.0 x 10^^ cm~^ are fixed by the SdH frequencies. For the 5-QL 
film (Fig. S5), where SdH oscillations were not detected, we assume rii — n2 — 'i y. 10^^ cm~^ following the trend in 
Us vs t for the surface channel [see Fig. 3(a)]. 

A three-band model analysis for the 8-QL film gives the mobilities /j,si = 1600 cm^V^^s^^ and /is2 = 1000 
cm^V~^s~^ for two different surfaces and /j,3 = 161 cm^V~^s~^ for the bulk carriers (Fig. S4). Note that this 
combination of carrier densities and mobilities for the three channels reproduce the curved Ryx{B) behavior as well 
as the B = T value of Rg, but the observed B dependence of Rs is stronger than the three-band model prediction; 
this is understandable, because the simple multi-band model employed here neglects the intrinsic magnetoresistance 
{i.e. the B dependence of the mobility) which is known to exist in Bi2Se3. 

In contrast, the Ryx data for the 5-QL film is almost perfectly linear in B, and to reproduce this linear behavior 
while keeping the consistency with the Rs(0) value, one needs to consider a very low mobility for the surface channel. 
For example, a combination of the surface mobility /ii^2 — 50 cm^V~^s^^ and the bulk mobility /is — 140 cm^V^^s^^ 
gives reasonable reproductions of Rs{0) and Ryx{B), as shown in Fig. S5. On the other hand, if one assumes a gradual 
change in the surface mobility across tc and hypothetically put fj,i^2 = 500 cm^V~^s~^, it is impossible to make Rs{0) 
to be consistent with the data nor to make Ryx{B) linear, as shown in Fig. S6. Therefore, one can conclude that the 
surface mobility indeed dropped by more than one order of magnitude when t was reduced from 8 to 5 QL. 



S5. Quantum confinement effects 



A. Formation of quantum-vifell states in ultrathin films 

In thin films, the effects of quantum confinement modify the 3D bulk energy states to 2D quantum-well (QW) 
states when the film thickness becomes comparable to the de Broglie wavelength of the carriers. In the case of very 
thin films, it is usually difficult to distinguish between a genuine SS and a 2D QW state if one simply measures the 
angular dependence of the SdH oscillations. Fortunately, we were able to measure SdH oscillations in all films with 
t > 6 QL, including those in which bulk states are certainly in the 3D limit. As an example. Fig. S7 shows the 2D 
character of SdH oscillations measured in a 120-QL film. It is worth noting that in none of the measured films we 
saw SdH oscillations coming from the 3D bulk states. This is most likely because of unfavorable conditions for their 
observations, that is, the relatively low mobility and a high bulk-carrier concentration together make the oscillation 
amplitude to become smafi 0. The systematic measurements of SdH oscillations in a wide range of thickness make 
it possible to reliably distinguish the SSs from bulk states or the 2D QW states by tracing them from thick samples. 



B. Surface band bending and resulting 2D electron gas 

The surface bend bending can create a 2D electron gas (2DEG) on the surface of TIs, and this could be responsible 
for the observed 2D SdH oscillations instead of the topological SSs. Spectroscopically, coexistence of the topological 
state and a 2DEG state on the surface of Bi2Se3 has been observed by ARPES a few hours after cleaving H, 5^ . 
which is believed to be due to gradual formation of defects (Se vacancies) on a freshly cleaved surface in ultrahigh 
vacuum. Figure S8 reproduces the energy dispersions of the SS and the 2DEG state measured along the K — T — K 
direction in Ref. [38| . Solid lines are our fittings of the equations shown in Fig. S8 to the experimental data. 

Now we show that it is not possible to ascribe the observed 2D SdH oscillations to the 2DEG state. The frequency 
of the SdH oscillations in the 10-QL film shown in Fig. 2 of the main text corresponds to kp — 5.7 x 10^ cm~^. As 
discussed in the main text, if the SdH oscillations are coming from the topological SS, this kp gives the estimate for 
Ep of 160 meV measured from the Dirac point. On the other hand, if this kp is associated with the 2DEG state 
shown in Fig. S8, Ep should be 82 meV from the bottom of this 2DEG state and the corresponding density of the 
2DEG, n2DEGi is 5.2 x 10^^ cm~^. One should remember that for this 2DEG state to form, the conduction band must 
be bent downward near the surface, which dictates that the Ep of the bulk state measured from the bulk conduction 
band edge must be smaller than 82 meV. This means that the maximum possible bulk carrier density is 5 x 10^^ 
cm^'^ if we use the effective mass m* = 0.13mo for the bulk conduction band [1^. (A larger bulk Ep would mean 
that the band is bent upward near the surface, which precludes the formation of the 2DEG state from the beginning.) 

However, the high-field slope of the Ryx{B) data, which simply gives a sum of the densities of electrons from all 
bands 0], points to the existence of a total of 4 x 10^^ cm^'^ of electrons in our 10-QL film; since the above 2DEG 
density corresponds to the volume density of only 5.2 x 10^^ cm^'^, the maximum possible bulk carrier density of 5 x 
10^* cm~'^ is not enough to account for the Hall data. [Note that the density of Dirac electrons is 1 x 10^'^ cm"^ at 
most in the above 2DEG scenario (see later), and hence the topological SS can only provide up to 2 x 10^^ cm^'^ of 
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carriers from the two surfaces, and this is still not enough.] Therefore, one cannot achieve a consistent understanding 
of the transport data in the 2DEG scenario. 

To illustrate the difficulty of the 2DEG scenario in a more quantitative manner, let us assume, for simplicity, the 
same band bending as observed in Ref. fsS] (i.e. the main 2DEG state located 40 meV below the bottom of the 
bulk conduction band and 90 meV above the conduction-band edge of the surface, which gives the total of 130 meV 
downward band bending, see the inset of the left panel of Fig. S9) [HI]. Then the position of the Dirac point (DP) 
is about 370 meV below the surface Fermi level, giving the concentration of Dirac electrons riDirac = 1 x 10^'^ cm^^. 
The bulk Ep should be 42 meV in this situation, and the corresponding bulk carrier density is ribuik — 1-8 x 10^^ 
cm~'^ ^]. Therefore, all carrier densities are fixed and the only parameters that can be tuned to reproduce i?s(0) and 
Ryx{B) are the mobilities of 2DEG, Dirac, and bulk electrons. Figure S9 shows the result of the best fitting, which 
actually does not fit the Ryx{B) behavior at all and obviously the slope of the calculated Ryx{B) is way too steep. 
This is essentially because the concentrations of both Dirac and bulk electrons are too low in the 2DEG scenario. As 
presented in the main text, the alternative scenario to ascribe the SdH oscillations to the topological SS and consider 
a slight upward band bending (which precludes the formation of the 2DEG state) naturally explains the relatively 
large bulk carrier density and gives a totally consistent description of our data. 

It is worth noting that the above-estimated Ep of 82 meV for the 2DEG state in the downward band-bending 
scenario is already significantly smaller than that seen in the ARPES data shown in Fig. S8, where Ep is about 200 
meV for the 2DEG state. This means that in our 10-QL film, even if one assumes that the 2DEG state were formed 
and were responsible for the SdH oscillations, the accumulated charge on the surface must be much smaller than in 
the case observed in ARPES experiments [sl], and hence the band bending, if existed, must be much weaker [s^ ]. 
In other words, the SdH data dictate that the band bending as strong as that shown in Fig. S9 cannot actually be 
conceived for our 10-QL film. 

In addition, we note that we did not observe any sign of progressive downward band bending in our transport 
data, even though our samples are not protected by any passivation layer and we applied a couple of different sample 
preparation procedures: some samples were transferred into an inert atmosphere immediately after the growth and 
were exposed to air for only several minutes before measurements, other samples spent hours in air during the contact 
preparation, and several samples have been remeasured months after their first measurement. So far no relaxations 
in transport properties have been detected in our MBE films, in contrast to what we had observed in cleaved bulk 
single-crystal samples fsl or what has been reported in the literature [JH^llB'l. It is possible that our films, which are 
of high-enough quality to show SdH oscillations, are inherently more stable than films of lower quality; namely, if the 
selenium-terminated surfaces have smaller number of defects and hence contain fewer dangling bonds, their chemical 
stability is expected to be higher [H^ and they are better protected against oxidation fsBI and the resulting charging 
effect. 



S6. Film quality across tc 

To confirm that the observed diminishment of the surface transport in ultrathin films is not caused by a lowering 
of the quality in thinner films, we have characterized the morphology of our films using atomic force microscopy. As 
shown in Fig. SIO, our ultrathin films remain essentially flat and smooth across tc (surface bumpiness of only ^^1 
QL). This observation suggests that it is not very likely that there is a drastic degradation in the film quality as a 
function of t. 

It is worth noting that in our films the sheet density of bulk carriers [shown in Fig. 3(a) of the main text by circles] 
changes smoothly but it is not proportional to t, which is similar to what was previously reported [23| : the observed 



trend suggests that the density of Se vacancies gradually increases in thinner films 23[ . 

S7. Weak antilocalization 

The fitting parameters of the WAL behavior measured in our films are shown in Table I. One can clearly see that 
both a and rapidly decrease below the critical film thickness. A similar tendency has also been observed in previous 
thickness-dependence studies although the decrease of a was less pronounced and more gradual, which is not 



a surprise for more metallic films where the suppression of the topological SS is less crucial. 

We also measured the angle dependence of the WAL behavior. Figure Sll shows the low- field magnetoresistance 
measured in the 10-QL film in tilted magnetic fields. The WAL response plotted as a function of BcosO clearly 
demonstrates its entirely 2D character and is well fit by the Hikami-Larkin-Nagaoka (HLN) formula [Eq. (1) of the 
main text]. 
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Thickness (QL) 




1/0 (nm) 


2 


0.007 


101 


3 


0.04 


153 


5 


0.5 


256 


8 


0.43 


409 


9 


0.45 


584 


1 n 

xu 




o 1 o 


15 


0.45 


573 


21 


0.5 


794 


32 


0.4 


405 


50 


0.57 


905 



TABLE I: Parameters of the WAL behavior described by the HLN formula. 
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FIG. SI: (a) Schematic picture of the Bi2Se3 crystal structure (two quintuple layers are shown), (b) Schematic picture of the 
van der Waals MBE growth on sapphire. 
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FIG. S2: Temperature dependence of the sheet resistance in a 2-QL film showing a diverging trend at low temperatures. Inset 
shows the temperature dependence of the Hall mobility /i — l/enRs- 




FIG. S3: SdH oscillations in dRxx/dB measured in films with various thicknesses. Insets show the Fourier transform of the 
data in the main panel. 
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FIG. S4: Fitting of a three-band model to Rs{B) (left panel) 
Symbols are experimental data. Solid lines are calculations. 
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and Ryx{B) (right panel) measured in the 8-QL film at 1.6 K. 




FIG. S5: Two-band model calculations of Rs{B) (left panel) and Ryx{B) (right panel) for the 5-QL film with the parameters 
shown in the left panel, together with the experimental data measured at 1.6 K. 
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FIG. S6: Two-band model calculations of Rs{B) (left panel) and Ryx(B) (right panel) for the 5-QL film with a hypothetical 
surface mobility of 500 cm'^V^^s^^ which is lower than that for t > tc but is still reasonably high. The experimental data 
measured at 1.6 K are shown by symbols for comparison. 
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FIG. S7: dRyx/dB (shifted for clarity) measured in tilted magnetic fields in the Bi2Se3 film with the thickness of 120 QL and 
plotted as a function of l/Bcos6, demonstrating the 2D character of oscillations with two frequencies. Inset shows the Fourier 
transform of the data at S = 0° indicating the two-frequency nature of the oscillations. 
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FIG. S9: Rs{B) and Ryx{B) measured in the 10-QL film (symbols) and their fitting (solid lines) using parameters constrained 
in a 2DEG scenario (see text). Inset in the left panel shows a schematic picture of the band bending considered for the 
calculations, where Ec is the bottom of the conduction band, Ev is the top of the valence band, Ep is the Fermi level, Eq is 
the main (first) 2DEG energy state appeared due to the confinement of the electron motion perpendicular to the surface of the 
film. 
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FIG. Sll: The WAL behavior measured for various magnetic field directions in the 10-QL film. 6 is the angle between the 
magnetic field and the surface normal. The dashed line is the fitting of the HLN formula. 



